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We show that inconclusive photon subtraction (IPS) on twin-beam produces non-Gaussian states 
that violate Bell's inequality in the phase-space. The violation is larger than for the twin-beam 
itself irrespective of the IPS quantum efficiency. The explicit expression of IPS map is given both 
for the density matrix and the Wigner function representations. 
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I. INTRODUCTION 

The twin-beam state (TWB) of two modes of radiation 
can be expressed in the photon number basis as 



oo 



|A))ah = v/l-A2 V A"|n,n),fc, 



E 

n=0 



(1) 



where A = tanh(r), r being the TWB squeezing parame- 
ter. TWB is described by a Gaussian Wigner function 

Wr{a,l3) = Aexp{-2A(|a|2 + |/3n 

+2B{aP + aP)}, (2) 

with A = A{r) = cosh(2r) and B = B{r) = sinh(2r). 
Since 10) is positive-definite, TWB are not suitable to 
test nonlocality through homodyne detection. Indeed, 
the Wigner function itself provides an explicit hidden 
variable model for homodyne measurements [ij, l2| ■ On 
the other hand, it has been shown [y] that TWB exhibits 
a nonlocal character for parity measurements. This is 
known as nonlocality in the phase-space since Bell in- 
equalities can be written in terms of the Wigner func- 
tion, which in turn describes correlations for the joint 
measurement of displaced parity operators. Overall, the 
positivity or the negativity of the Wigner function has a 
rather weak relation to the locality or the nonlocality of 
quantum correlations. 

In Ref. 1^ we have suggested a conditional measure- 
ment scheme on TWB leading to a non Gaussian entan- 
gled mixed state, which improves fidelity in the teleporta- 
tion of coherent states. This process, called inconclusive 
photon subtraction (IPS), is based on mixing each mode 
of the TWB with the vacuum in a unbalanced beam split- 
ter and then performing inconclusive photodetection on 
both modes, i.e. revealing the reflected beams without 
discriminating the number of the detected photons. 

A single mode version of the IPS, mapping squeezed 
light onto non-Gaussian states, has been recently realized 



experimentally jj|. Moreover, IPS has been suggested as 
a feasible method to modify TWB and test nonlocality 
using homodyne detection [2|. 

In this paper we address IPS as a degaussification map 
for TWB, give its explicit expression for the density ma- 
trix and the Wigner function, and investigate the nonlo- 
cality of the resulting state in the phase-space. 

The paper is structured as follows. In Section UTI we 
review nonlocality in the phase-space, i.e. Wigner func- 
tion Bell's inequality based on measuring the displaced 
parity operator on two modes of radiation. In Section 
mil we illustrate the IPS process as a degaussification 
map and calculate the Wigner function of the IPS state. 
The nonlocality of the IPS state in the phase-space is 
then analyzed in Section IIVI whereas in Section El we 
discuss nonlocality using homodyne detection, extending 
the analysis of Ref. [3. Section IVll closes the paper with 
some concluding remarks. 



II. NONLOCALITY IN THE PHASE-SPACE 

The displaced parity operator on two modes is defined 
as 

U{a,P) = Da{a)(-\f''Dl{a) 

®D,{[i){-lf'Dl{(3), (3) 

where a, /? e C, a and b are mode operators and Da{a) — 
expjaa''^ — aa] and Db{f5) are single- mode displacement 
operators. Parity is a dichotomic variable and thus can 
be used to establish Bell- like inequalities jaj. Since the 
two-mode Wigner function W{a, (3) can be expressed as 



W{a,l3) 



■n(a,/?). 



(4) 
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n(a, /?) being the expectation value of 11(0;, /3), the viola- 
tion of these inequalities is also known as nonlocality in 
the phase-space. The quantity involved in such inequali- 
ties can be written as follows 

B = n(ai,/?i) + n(a2,/3i) 

+n(ai,/?2)-n(a2,/32), (5) 

which, for local theories, satisfies the condition \B\ < 2. 




FIG. 1: Scheme of the IPS process. 

Following Ref. [l|, one can choose a particular set of 
displaced parity operators, arriving at the following com- 
bination 



B{J) = n(0, 0) + n(\/j, 0) 

+n(o, -Vj) - n{Vj, -Vl) , 



(6) 



which only depends on the positive parameter J, charac- 
terizing the magnitude of the displacement. If we eval- 
uate the quantity © in the case of the TWB, we find 
that it exceeds the upper bound imposed by local the- 
ories for a certain region of values {J,r), its maximum 
being S« 2.19 Q. 

On the other hand, the choice of the parameters lead- 
ing to Eq. (jnj is not the best one, and the violation of 
the inequality |S| < 2 can be enhanced using a different 
parameterization [6j. A better result is achieved for 

c{j) = n(\/j,-\/j) + n(-3\/j,-\/j) 

+n(\/j, 3Vj) - n(-3\/j, 3Vj) , (7) 

which, for the TWB, gives a maximum C sa 2.32, greater 
than the value 2.19 obtained in Ref. [l|. 

In the following Sections we will see that the violation 
of the inequalities |i3(J)| < 2 and |C(J)| < 2 can be 
enhanced by degaussification of the TWB. 



III. THE DEGAUSSIFICATION PROCESS 



The degaussification of a TWB can be achieved by sub- 
tracting photons from both modes [3jSl3- In Ref. |3| we 
referred to this process as to inconclusive photon subtrac- 
tion and showed that the resulting state, the IPS state 
ftps, can be used to enhance the teleportation fidelity 
of coherent states for a wide range of the experimental 
parameters. 

The IPS scheme is sketched in Fig. ^ The two modes, 
a and 6, of the TWB are mixed with the vacuum (modes 
c and d, respectively) at two unbalanced beam splitters 
(BS) with equal transmissivity r = cos^ (j>; the modes 
c and d are then revealed by avalanche photodetectors 
(APD) with equal efficiency rj. APD's can only discrimi- 
nate the presence of radiation from the vacuum. The pos- 
itive operator- valued measure (POVM) {no(?7). Hi (77)} of 
each detector is given by 



no(r;) = 5^(l-r/F|j)01, ni(77)=I-no(77), (8) 

3=0 



rj being the quantum efficiency. Overall, the conditional 
measurement on the modes c and d, is described by the 
POVM (we are assuming the same quantum efficiency for 
both photodetectors) 



noo(r?) = no,c(?7)®no,d(??) 

nio(?7) = ^iAv)'»^o.d{v) 
nii(?7) = ni,c(?7)«)ni,d(77) 



(9) 
(10) 

(11) 

(12) 



When the two photodetectors jointly click, the condi- 
tioned output state of modes a and b is given by 



J 



£{R) = 



PiiirAiV) 



Ti-cd[Uac{(l)) ® Um{(J))R® |0)cc(0| ® |0)drf(0| UlM) ® UlM la ® lb ® Ilii{v)] , (13) 



where Uac{4') = exp{— 0(a^c — ac^)} and Ubd{4') ^-re the evolution operators of the beam splitters and R the density 
operator of the two-mode state entering the beam splitters (in our case R ~ Qtwb ~ \-^))abba{{M)- The partial trace 
on modes c and d can be explicitly evaluated, thus arriving at the Kraus decomposition of the IPS map. We have 



£iR) = 



1 



Pii{r,(j),T]) 



J2 ^^pi^: V) Mpgi<P) R A'4(0) m,(0, 7j) 



p,q=l 



(14) 



with 



mpicj),!]) = 



t&n^Pcj) [1- {l-rj)P] 



and 



(15) 



Mpq{(j))= aP6'?(cos(/))'^'"+^^^ 



(16) 



and 



Pii{r,(j>,'n) = Trab[£{R)] 



(17) 



is the probability of a click in both detectors. 

Now, in order to investigate the nonlocality of the state 



J 



Pips ~ f (ptwb) in the phase-space, we explicitly calculate 
its Wigner function, which, as one may expect, is no 
longer Gaussian and positive-definite. 

The state entering the two beam splitters is described 
by the Wigner function 



W^-\a,f3,C,O^Wr{a,(3)—cxp{-2\C\^-m^} , 



(18) 



where the second factor at the rhs represents the two vacuum states of modes c and d. The action of the beam 
splitters on W^'"' can be summarized by the following change of variables 



a —^ a cos 4> + C sin (j) , C ^ C cos cj) — a sin , 
(3 ^ P cos -t- ^ sin , (, ^ £, cos — /3 sin (/) , 



(19) 
(20) 



and the output state, after the beam splitters, is then given by 

4 



w^;7'("'/3'C-o 



where 



and 



exp {-a|Cp + {v + 2B£,sm'^ 0)C + {v + 2BCsin2 0)C} 
f \(3\^) + 2B cos^ (j) (aP + a^)} 



Wr,4,ia,(3) = -^ exp{-6(| 



a = a(r, (/)) = 2(Asin^ (/) -I- cos^ (/)), 

h = h{r, (f)) = 2{A cos^ (j) + sin^ </>) , 

V = w(r, 0) = 2cos0 sin0 [(1 — A)a 4- _B/3], 

w = w{r,(j)) =2 cos (j) SUV (f)[[l — A) P + B a] 

I 



(21) 



(22) 



(23) 
(24) 
(25) 
(26) 



At this stage conditional on/off detection is performed 
on modes c and d (see Fig. ^|. We are interested in 
the situation when both the detectors click. The Wigner 
function of the double click element 1111(77) of the POVM 
(see Eq. (O) is given by [1I3 



t^,(C,0 ^ w^[nii(,7)](c,6 

= ^{i-Q.(C)-Q.(0 

+Qn{OQn{0}^ 



with 



'>,(^) 



2-77 



exp • 



_2rj_ 
2-77 



(27) 



(28) 



Using Eq. (|13|1 and the phase-space expression of trace 
(for each mode) 

Tr[0i02] = TT / d^zW[0i]{z)W[02]{z) , (29) 

Jc 
Oi and O2 being two operators and Ty[Oi](z) and 

W^[02]('Z) their Wigner functions, respectively, the 

Wigner function of the output state, conditioned to the 

double click event, is then given by 



Wr,4,,r,ia, P) 



pii(r,0,77) 



(30) 



where 



fr,4,,r,{a, P) =7r^ 



d\d'^^WrAc^,(3)Y,^Gi%^{a,P,<:,0 



(31) 



TABLE I: 



la a 1 

2 a+TT— a -T^- 

' 2 — 77 2 — 77 

2 — 77 2 — 77 ^2 — ri^ 

and pii(r, 0, if) is the double-click probability (|17fl . which can be written as function of fr,(j),ri{oL^ (3) as follows 

Pn(r,</.,77)=^2 /■ d^ad^l3fr,4>,v^a,P)- (32) 

The quantity G^-*], (a, /3, C, C) appearing in Eq. (|^ is 

GS,r,("'AC,0 = exp{-x,|C|' + (i' + 2B^sin2 0)C + (iJ+2B?sin2 0)C} 



xexp{-y,|e|2+< + w^} , (33) 
and the expressions of Cj{rf), Xj = Xj{r, (p, rj) and yj = yj{r, </>, 77) are given in Tabic p] 

The mixing with the vacuum in a beam sphtter with foUowed by an ideal {i.e. efficiency equal to 1) on/off 

transmissivity r followed by on/off detection with quan- detection. Therefore, the state H3l)|l can be studied for 

turn efficiency n is equivalent to mixing with an effective 77 = 1 and replacing r with Teff. Thanks to this substitu- 

transmissivity [3J tion, after the integrations we have 

Teff = Tcff (0, 77) = 1 - 77(1 - t) (34) 



1 4 

/r,0,,)(a,/3) = — y" 



16Q(77) 



X exp{-(fe - /j)|ap - (6 - .gj)l/3|' + {^Bt.s + hj){al3 + aP)} (35) 



and 



pii{r, (/>, 77 = > ., -^ " w, X 7:7^7 T^— > (36) 



where we defined 



fj = f,{r,cj),i^) = N,[x,{l-Af + W\l-A){l~T,s) + V,B^], (37) 

g, = g^ir, 0, 77) - iV, [x^B^ + 4B\1 - A)(1 - r,ff) + 77,(1 - Af] , (38) 

h, = h,{r, cj,, Tj) = N, [{x, + y,)B{l - A) + AB{B^ + (1 - Af){l - Tes)] , (39) 

In this way, the Wigncr function of the IPS state can be rewritten as 

w ! m w ! ,, ^ 4Q(77)J^g.^(a,/3) 



where we introduced 



K. 



U) 



';ja,(3) = exp{/,|a|^ +5,|/3|^ + /i,(a/3 + «/3)} 



The state given in Eq. (|1T|| is no longer a Gaussian state. 



(42) 
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FIG. 2: Plot of B{J) given in Eq. ® for J = 10'^ The 
dashed line is i3r^™^'(J), while the solid lines are B!^^ hj) 
for different values of Tes (see the text): from top to bottom 
Tcff = 0.999, 0.99 and 0.9. When TcS = .999, the maximum 
of S^'^^'(J) is 2.23. The lower plot is a magnification of the 
region < r < 0.11 of the upper one. Notice that for small r 
there is always a region where i3r^^^'(J) < B'^^^hJ). 



IV. NONLOCALITY OF THE IPS STATE 
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In this Section we investigate the nonlocality of the 
state H41() in phase-space using the quantity B given in 
Eq. ((nj, referring to both the parameterizations B{J) 
(see Eq. (jHJ) and C{J) (see Eq. 0). 

As for a TWB, the violation of the Bell's inequality is 
observed for small r []|. From now on, we will refer to 
B{J) as B),^*^'(J) when it is evaluated for a TWB ©, 
and as B'^^JUj) when we consider the IPS state (|^ . 

We plot S7"=^)(J) and S^|^^^(J) in the Figs. HandElfor 
different values of the effective transmissivity Tcs and of 
the parameter J: for not too big values of the squeezing 
parameter r, one has that 2 < B'~^J^^\j) < B'~^^^^^^{J) . 
Moreover, when TeS approaches unit, i.e. when at most 
one photon is subtracted from each nrode, the maxinrum 
of B]^^^l is always greater than the one obtained using a 
TWB. A numerical analysis shows that in the limit t^s -^ 
1 the maxinrum is 2.27, that is greater than the value 2.19 



FIG. 3: Plots of B[,J) given in Eq. @ as a function of 
the squeezing parameter r for different values of J: (a) J = 
5 10"^ (b) J = 10"^ and (c) J = 2 10"\ In all the plots the 
dashed line is ^^"'''(J'), while the solid lines are -B<'^^^(J) 
for different values of TcS (see the text): from top to bottom 
Tcff = 0.999,0.9,0.8,0.7 and 0.5. Notice that there is always 
a region for small r where ^^"^'(J) < B^'^^^l{J). When 

Tcff 

theoneofB<™'"(J) 



0.999 the maximum of B^^^^'(J) is always greater than 



obtained for a TWB [l|. The Hmit t^s —^ 1 corresponds 
to the case of one single photon subtracted from each 
mode 0jl3- Notice that increasing J reduces the interval 
of the values of r for which one has the violation. For 
large r the best result is thus obtained with the TWB 





FIG. 4: Plots of C{J) given in Eq. as a function of the 
squeezing parameter r for J = 1.6 10~*. In all the plots the 
dashed line is Cr"'"(J), while the solid lines are C^]^^^(J) 
for different values of t^s (see the text): from top to bottom 
Tcff = 0.999, 0.99, 0.95, 0.9 and 0.8. When t^s = 0.999 the 
maximum 



ofQX^(J)is2.40. 



since, as the energy grows, more photons are subtracted 
from the initial state Q. Since the relevant parameter 
for violation of Bell inequalities is Tos , we have, from Eq. 
(|34|l . that the IPS state is nonlocal also for low quantum 
efficiency of the IPS detector. 

The same conclusions holds when we consider the pa- 
rameterization of Eq. (0. In Fig. E|we plot G^^^'^J) 
and C<|^^^(J), i.e. C{J) evaluated for the TWB and the 
IPS state, respectively. The behavior is similar to that of 
B{J), the maximum violation being now C^^^^^( J) = 2.40 
for Teff == 0.999 and J = 1.6 10""'. 

Finally, notice that the maximum violation using IPS 
states is achieved (for both parameterizations) when Tcfr 
approaches unit and for values of r smaller than for 
TWB. 



FIG. 5: Plots of S given in Eq. 14311 as a function of tanh(r-) 
for different values of T^a and for ideal homodyne detection 
(i.e. with quantum efficiency ryn = 1): from top to bottom 
Tcff = 0.99, 0.95, 0.90, 0.80 and 0.70. 




FIG. 6: Plots of S given in Eq. 14311 as a function of 
tanh(r) with t^s = 0.99 and for different values of the ho- 
modyne detection efficiency t^h: from top to bottom t^h = 
1, 0.95, 0.90, 0.85 and 0.80. The maximum of the violation de- 
creases and shifts toward higher values of r as 77H decreases. 
For smaller values of Teff the violation is furtherly reduced. 



V. NONLOCALITY AND HOMODYNE 
DETECTION 

The Wigner function Wr,0,,,(a,/3) given in Eq. l|^ 
is not positive-definite and thus pips can be used to test 
the violation of Bell's inequalities by means of homodyne 
detection, i.e. measuring the quadratures X{) and x^p of 
the two IPS modes a and 6, respectively, as proposed in 
Ref. [2|. In this case, if one discretizes the measured 
quadratures assuming as outcome +1 when a; > 0, and 
— 1 otherwise, one obtains the following Bell parameter 

+E{d2,^i)-E(§2,V2), (43) 

where "dj and ipj are the phases of the two homodyne 
measurements at the modes a and h, respectively, and 



Ei'&j.Vk) 



dx^j dx^^ sign[a;^^. x^J P{xi)- , x^ J , 



(44) 



P{x^-,Xip^) being the joint probability of obtaining the 
two outcomes x^. and x^^^ [3. As usual, violation of 
Bell's inequality is achieved when \S\ > 2. 

In Fig. El we plot S for i?i = 0, i?2 = ^f^^^i = 
— 7r/4 and (p2 = 7r/4: as pointed out in Ref. [21, the 
Bell's inequality is violated for a suitable choice of the 
squeezing parameter r. Notice that when Teff decreases 
the maximum of violation shifts toward higher values of 
r. 

As one expects, taking into account the efficiency % 
of the homodyne detection furtherly reduces the viola- 
tion (see Fig. |SJ). Notice that, when ?/h < 1, violation 
occurs for higher values of r, although its maximum is 
actually reduced: in order to have a significative viola- 
tion one needs a homodyne efficiency greater than 80% 
(when Teff = 0.99). 



VI. CONCLUDING REMARKS 

In this paper we have shown that IPS can be used 
to produce non-Gaussian two-mode states starting from 



a TWB. We have studied the nonlocality of IPS states 
in phase-space using the Wigner function. As for the 
improvement of IPS assisted teleportation Q, we have 
found that the nonlocal correlations are enhanced for 
small energies of the TWB (small squeezing parameter 
r). Moreover, nonlocality of pips is larger than that of 
TWB irrespective of IPS quantum efficiency. 

Since the Wigner function of the IPS state is not pos- 
itive definite, we have also analyzed its nonlocality using 
homodyne detection. In this case violation of Bell's in- 
equality is much less than in the phase-space, and is fur- 
therly reduced for non unit homodyne efficiency % < 1. 
However, this setup (IPS with homodyning) is of partic- 
ular interest, since it can be realized with current tech- 
nology achieving a loophole-free test of Bell's inequality 



Q 



On the other hand, the experimental verification of 
phase-space nonlocality is challenging, due to the difficul- 
ties of measuring the parity, either directly or through the 
measurement of the photon distribution. On the other 
hand, the recent experimental generation of IPS states 
[j| is indeed a step toward its implementation. 
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